Resonant activation in the presence of non-equilibrated baths. 
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We study the generic problem of the escape of a classical particle over a fluctuating barrier under the influence 
of non-Gaussian noise mimicking the effects of not-fuUy equilibrated bath. Our attention focuses on the effect 
of the stable noises on the mean escape time and on the phenomenon of resonant activation (RA). Possible 
physical interpretation of the occurrence of Levy noises in the system of interest is discussed and the connection 
between the Tsallis statistics and the Fractional Fokker Planck Equation is addressed. 
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I. INTRODUCTION 

In a classical Langevin approach the influence of the bath of 
surrounding molecules on a Brownian particle is described in 
terms of a mean-field, time-dependent stochastic force which 
is commonly assumed to be white Gaussian noise. That pos- 
tulate is compatible with the assumption of a short correla- 
tion time of fluctuations, much shorter than the time-scale of 
the macroscopic motion and assumes that weak interactions 
with the bath lead to independent random variations of the pa- 
rameter describing the motion. In more formal, mathematical 
terms Gaussianity of the state-variable fluctuations is a con- 
sequence of the Central Limit Theorem which states that nor- 
malized sum of independent and identically distributed ran- 
dom variables with finite variance converges to the Gaussian 
probability distribution. If, however, after random collisions 
jump lengths are ruled by broad distributions leading to the 
divergence of the second moment, the statistics of the process 
changes significantly. The existence of the limiting distribu- 
tion is then guaranteed by the generalized Levy-Gnedenko 
limit theorem. According to the latter, normalized sums of 
independent, identically distributed ii.i.d) random variables 
with infinite variance converge in distribution to the Levy 
statistics. At the level of the Langevin equation. Levy noises 
are generalization of the Brownian motion and describe re- 
sults of strong collisions between the test particle and the 
surrounding environment. In this sense, they lead to differ- 
ent models of the bath that go beyond a standard "close-to- 
equilibrium" Gaussian description. 

As documented elsewhere |2, 3], not fully thermalized sys- 
tems or systems driven away from the equilibrium can man- 
ifest interesting physical properties. In particular, such sys- 
tems may exhibit large energy fluctuations with probabilities 
higher than those predicted by the Gaussian statistics. Non- 
equilibrated heat reservoir can be thus considered as a source 
of non-Gaussian noises. Formalisms that give physical back- 
ground for this phenomena are based on the idea of nonexten- 
sive thermodynamics, established on a Tsallis statistics 0], 
and the Fractional Fokker Planck Equation (FFPE) \3\. The 
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latter (FFPE) have been successfully applied 

jisain for 

describing anomalous diffusion processes, for which the non- 
local character brought about by the stable Levy noise leads to 
the replacement of the local spatial derivatives in the diffusion 
term of the Fokker-Planck equation by a fractional derivative. 

In contrast to the fractional calculus, the Tsallis statistics 
approach to stable. Levy-type fluctuations is based on the 
introduction of a special entropy form ^ 1IQ1I from which 
examined thermodynamic al quantities can be derived 
The main property of the Tsallis entropy is its nonextensivity, 
i.e. the entropy of the system containing two noninteracting 
subsystems is different from the sum of the subsystems' en- 
tropies typical for equilibrium Gibssian ensembles and Gaus- 
sian measures. 

Various examples of everyday life phenomena and a wide 
extent of experimental observations f^, |3 fl3ll show exis- 
tence of long-range coiTelations, disorder, cooperativity and 
deviations from the Gaussian statistics, thus suggesting a 
strong need of study more general probability distribution 
than just Gaussian ones. In particular, a-stable distributions 
have been observed in anomalous dynamics and strange ki- 
netics in amorphous semiconductors and glassy systems |8j]. 
Levy-flight models turn out to be adequate for the description 
of transport in heterogeneous catalysis, self-diffusion in mi- 
celle systems and analysis of geophysical data L14. .15^ .16il . 
Related models have been also applied in financial modeling 
and analysis of economic time-series fvfJT^. Among various 
aspects of the Levy-type variables and processes, a problem of 
special interest is numerical generation of stable distributions 
and simulation of a-stable integrals and stochastic differential 
equations |19, 20, 21, 22, 23, 24]. As broadly discussed in 
literature 1 20, 22, 24], crucial difficulties in solving stochastic 
differential equations with stable measures are caused by the 
noisy term, which allows for larger fluctuations with higher 
probabilities than Gaussian distributions. Numerical meth- 
ods 1 17, 20] for such equations are more sophisticated than for 
differential equations |25] and for stochastic differential equa- 
tions with Gaussian noises |26]. In particular, the nonexis- 
tence of variance for stable variables makes the problem much 
more complicated to tackle both, numerically II7II20I1 and an- 
alytically 127*, "28"]. It turns out, however, that with the use 
of suitable statistical estimation techniques, computer simu- 
lation procedures and numerical discretization methods it is 
possible to construct relevant approximations of stochastic in- 
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tegrals with stable measures as integrators f20l l2ir . 

In this communication we study statistical properties of the 
generic system describing passage of a classical particle over 
a fluctuating potential barrier. The system is coupled to a non- 
Gaussian bath modeled by the Levy stable noise. We present 
numerical results for the mean first passage time of the par- 
ticle over the barrier, assuming a linear potential subject to 
Markovian dichotomous fluctuations. Our model belongs to 
the class of "on-off" models discussed in a paper by Doer- 
ing and Gadoua |29] and further analyzed by Boguna et al. 
1 30]. A distinctive characteristics of this model is that part 
of the time the barrier is either switched off (i.e. it becomes 
flat), or the switching is performed between the barrier and a 
well, so that the particle can essentially roll rather than climb 
during these times. The main difference between the model 
considered here is a form of driving, additive fluctuations. 
Whereas in the previous papers mainly white Gaussian noises 
have been considered |29, 31, 32, 33, 34], here additive noises 
are assumed to be a-stable ]21], that might arise from the con- 
tact with not fuUy equilibrated bath. 

The problem of resonant activation (RA) f2^ examined in 
this paper is an example of processes manifesting constructive 
role of noise jHI. The article deals with a modified version 
of the model proposed by Doering and Gadoua f29']. Section 
II presents the model and poses the problem to be studied. In 
Sections III and IV general considerations of Levy noises in 
physical systems are addressed and the problem of underlying 
thermodynamic interpretation is discussed. Results of sim- 
ulations with a short note on numerical methods applied for 
generating stable variables and integration of stochastic differ- 
ential equations with stable measures are included in Section 
V. The paper is closed with the concluding remarks. 




FIG. 1: A model potential studied in the paper. The barrier height 
fluctuates dichotomously between the values H±. A particle starts 
its diffusive motion at a reflecting boundary a; = and continues 
until the absorption at x = 1. 



rate 7 (cf. Fig.Q 



V±{x) = H±x, g = 



(2) 



Both C and rj noises are assumed to be statistically indepen- 
dent. 

The initial condition for Eq. Q is 



a;(0) = 0, 



(3) 



II. GENERIC MODEL SYSTEM 

We consider an overdamped Brownian particle moving in 
a potential field between absorbing {x = 1) and reflecting 
(x = 0) boundaries (cf. Fig.[3, in the presence of noise that 
modulates the barrier height. 

Time evolution of a state variable x{t) is described in terms 
of the Langevin equation 



dx 
'dt 



= -v'ix)+g7m + at) 

= -Vi{x)+C{t), 



(1) 



where prime means differentiation over x, ({t) is a white 
Levy process originating from the contact with not-fuUy 
equilibrated bath and ri{t) stands for a Markovian dichoto- 
mous noise of intensity g taking one of two possible val- 
ues ±1. Autocorrelation of the dichotomous noise is set to 
iivit) - {V))ivit') - iv))) = exp(-27|i - t'\). For sim- 
plicity, throughout the paper a particle mass, a friction coef- 
ficient and the Boltzman constant are all set to 1 . The time- 
dependent potential V±{x) is assumed to be linear with the 
barrier switching between two configurations with an average 



i.e. initially particle is located at the reflecting boundary with 
equal choices of finding a potential barrier in (±) configura- 
tions 



P{H+,t^ 0) = P{H_,t^ 0) 



1 



(4) 



The quantity of interest is the mean first passage time (MFPT) 



MFPT= dx [p-{x,t) +p+{x,t)]dt, (5) 
Jo Jo 

i.e., the average time which particle spends in the system be- 
fore it becomes absorbed. Within the proposed approach, the 
MFPT is estimated as a first moment of the distribution of first 
passage times (FPT) obtained from the ensemble of simulated 
realizations of the stochastic process in question. Otherwise, 
for Levy flights described by Eq. Q, the MFPT may be calcu- 
lated after solving a relevant deterministic fractional Fokker- 
Planck equation ] 2fi] (FFPE) for the distribution function 



dp±{x,t) 
dt 



d_ 

dx 



dV±{x) 
dx 



p±{x,t) 



+ DV°'p±{x,t) 



(6) 
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In the above FFPE p±{x, t) are probability density functions 
(PDFs) for finding a particle at time t in the vicinity of x, 
while potential takes the value V±{x). The fractional deriva- 
tives V" is understood in the sense of the Fourier Trans- 
form L37.1 



J Svr 



(7) 



with a = 2 corresponding to the standard Brownian diffu- 
sion. The coefficient D denotes the generalized diffusion co- 
efficient with the dimension |36] [D] ~ cm^sec^^ and can 
be related to the parameter a characterizing the width of the 
PDF (see below). 



m. LEVY TYPE VARIABLES 

The a-stable variables are random variables for which the 
the sum of random variables is distributed according to the 
same distribution as each variable, i.e. 



aXi + bX2 = cX + d, 



(8) 



analytical expressions corresponding to Eq. (|9} can be given 
only in few cases. In particular, for a — 2, (3 — 0, C, is a 
Gaussian variable with the probability density 

1 f 

L,,o(x;a,M) = ^exp(^ 4^ j ' ^^2) 

whereas a = 1, /3 = and a — ^, /5 = 1 yield Cauchy- 

1 



Lifi{x]a,y) = - 



TT (x — flY + (T^ 



and Levy-Smirnoff {x > y)- 



X exp 



2(.T - ^I) 



(13) 



(14) 



distributions, respectively. 

Generally, for (3 — ^ — PDFs are asymmetric and for 
/3 = ±landae (0,1) they are totally skewed ilUl . 



where = denotes equality in a distribution sense. Real con- 
stants c, d in Eq. ^ allow for rescaling and shifting of the 
initial probability distribution. 

The characteristic function of the probability distribution 
that fulfills ^ can be parameterized in various ways |24]. In 
the usually chosen La^p{C,\(j, ji) 21, 22] parameteriza- 
tion, a characteristic function of the Levy type variables is 
given by 



0(fc) ~ exp 



-a"|fc|" (l - z/?sign(/c) tan ^) 



i(3kcr°' tan 



na 



, for a 1, 



m 



exp 



-a\k\ 1 + i/3-sign(fc)ln|A:| 



for a = 1, 



(9) 



with a e (0, 2], (3 £ [-1, 1], a € {Q,oo), ^ G (-c»,c») 
and (p{k) defined in Fourier space 



(10) 



Parameter a is called the stability index, f3 describes skew- 
ness of the distribution, a is responsible for its scaling of and 
/i is a location parameter The above parameterization (|9j is 
continuous, in the sense that 



lim a" (|/fc|"sign(fc) - fc)tan 

a— >1 2 



—(T\k\sign{k) In \k\ 

TT 

(11) 

for every a and k. 

Although the PDFs for stable variables L{() are known to 
have an asymptotic power-law behavior L{() ~ ICI"*'"^^''' 



IV. LEVY NOISES AND TSALLIS THERMODYNAMICS 

Classical equilibrium thermodynamics can be derived by 
use of the variational principle for the Boltzman-Gibbs- 
Shannon entropy issil 



S[p{x)] 



p{x) hi[ap(x)]dx, 



(15) 



where p{x) is a probability distribution of a given quantity. 

For example, x may be associated with a single jump length 
of a free Brownian particle fTol l3^ so that p{x) stands for 
a single jump length distribution. Optimizing S[p{x)] with 
additional constraints imposed on p{x) 



and 



p{x)dx = 1 



x'^p{x)dx = 



(16) 



(17) 



leads to a Gaussian distribution vix) used in the description 
of a normal Brownian diffusion JSj] 



p{x) exp(-fea;^)/Z, 



(18) 



Here Z — exp{—bx^) = y^n/b is a partition function 
and 6 is a Lagrange multiplier, b = 1/(2(t^). 

As an extension of the above, a generaUzed g-entropy sug- 
gested by Tsallis |4i] 



Sg\p{x) 



i-/_°°^M^r 

9-1 



dx 



(19) 
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with a supplementary condition 

/oo 
x^[ap{x)Yid{x/a) =a'^ (20) 
-OO 

leads to O 

p,(:r) = [l-fe(l-g)x2]i/(i-9)/Z„ (21) 

with Zq ^ /^[l - fe(l - g)a;2]i/(i-9)(fa;. In the hmit of 
q — > reproduces Eq. ( I18> and the norm-constraint can 

be satisfied only for q < 3. A^-folded distribution Pq{x^ N) 
approximating probability distribution for macroscopic diffu- 
sion is determined by the general central limit theorem. Con- 
sequently, for q < 5/3, the asymptotic {N — > oo) distribution 
is Gaussian whereas for q > 5/3 it approaches the Levy statis- 
tics iHIH 

Pqix, N) « i„((\/6/iVi/")x) (22) 
with the stability index a given by 

a = ^— I (5/3 < g < 3). (23) 

The variational principle can be also applied to the energy dis- 
tribution i40ll which subject to the generalized condition 

^[p(£0]'£. = (24) 

i 

results in 

Pq{e.,)^[l-b{l-q)e,f^'"'\ (25) 

Alternatively, expression Eq. (I25> can be derived by assuming 
subordination of the Gibssian exponential energy statistics to 
the statistics of fluctuations in temperature T (and hence in the 
parameter 6 = 1/T) |41] 

/•oo 

Pq{ei) cx / exp{-be)f{b)db (26) 
Jo 

In particular, a gamma-distributed temperature leads to a Levy 
distribution Pq{e) = La{be) with the stability index given by 
the relative variance of fluctuations in & |41]. This observation 
shows a possible link of the Levy distributions to models of 
nonequilibrated baths. 

V. PHYSICAL INTERPRETATION 

By definition, in the case of static barrier height and with 
the noise ( uncorrected at different times and obeying the 
Levy statistics, the overdamped Langevin equation Eq. Q de- 
scribes Levy flights 06[l37ll in a constant force field. For long 
times, the trajectory x{t) behaves as 

x{t) KiHot+ [ dsCis) (27) 
Jo 



and yields the Levy stable distribution f24','37l in the posi- 
tion of the particle. In consequence, if the first moment exists, 
i.e. for 1 < a < 2, mean value of x{t) grows linearly with 
time, {x{t)) — Hot whereas the mean square displacement 
becomes {[x{t) - {x{t))]'^) = 2D ^ 2a^ only for a = 2 
when the finite second moment of £(C) exists. Thus the gen- 
eralized Einstein relation connecting the first moment in the 
presence of constant force Hq to the second moment in the 
absence of force {x{t))Ho = Ho{x^{t))o is recovered only in 
the Brownian limit a — 2 with the noise amplitude a related 
to the diffusion coefficient a — D^l^ . Obviously, for noises 
with diverging mean square displacement, the classical fluctu- 
ation dissipation theorem is violated and the Einstein relation 
does not longer hold QHlil. 

When analyzed from the perspective of the continuous time 
random walks (CTRW), Levy flights characterize walks with 
a Poisson waiting time and a Levy distribution of the jump 
length. The scaling nature of the jump length PDF leads then 
to a clustering of the Levy flights visible via interruption of 
the local motion by occasional long sojourns on all length 
scales Is!]. This fractal character of the Levy flight- trajectory 
can be contrasted with subdiffusive CTRW |5, 9] that fills the 
two-dimensional space completely and features no clusters. 
In such a case, however, the time intervals between consec- 
utive steps are governed by the power-law waiting time dis- 
tributions that lead to a sublinear dependence (x(t)) w t" 
with V denoting the power-law index of the waiting time PDF. 
Therefore, subdiffusive CTRW in a constant external force 
acting along the x direction perfectly satisfies the fluctuation- 
dissipation theorem |5, 9] which holds also in the correspond- 
ing fractional Fokker-Planck equation framework. 

Diverging mean-square displacement cannot be valid for a 
particle with non-diverging mass. In fact, for massive parti- 
cles, a finite velocity of propagation exists making very long 
instantaneous jumps impossible. For that reason, the dilemma 
of diverging mean-square displacement in the Levy flight can 
be overcome by CTRW version of Levy walks with a suitable 
time cost penalizing long jumps. Nevertheless, in many phys- 
ical systems of interest diffusion of state variable x(t) with 
diverging second moment does not violate physical principles 
and is a legitimate way of physical modehng LI 4.. 15. ,43. .44.1 . 



VI. STOCHASTIC DIFFERENTIAL EQUATIONS WITH 
STABLE NOISES 

Instead of solving (|6j the corresponding Langevin equa- 
tion ([0 can be simulated by use of the appropriate numerical 
methods. Position of the particle is then obtained by direct 
integration of Eq. Q 

x{t) = - f [V'{x{s))-grj{s)]ds 

Jto 

+ / dLa.pis) 

Jto 

= - f V±ixis))ds+ f dL^^p{s). (28) 

Jto Jto 



5 



In general the Lo..f3 measure in Eq. i28t can be 

approximated by 



f f{s)dLo.As) « f{^^sWo.,P iV^s, {i + l)As)) 
^ ^ /(zAs)Asi/^„ (29) 

i=0 

where is distributed with the PDF La^p{<;] cr = 1, /i = 0), 
N/S.S = t — to and AfQ_;5([iAs, (i + l)As)) is the measure of 
the interval [iAs, {i + l)As). 

Random variables ^ corresponding to the characteristic 
functi on (l9 t can be generated using the Janicki-Weron algo- 
rithm ll7ll20lE2il . For a ^ 1 their representation is 



sin(a(y + Ca^p)) 



(cos(y))i 



W 



with constants B,C,D given by 



(30) 



Ba,i3,a.,i = M - ^cr" tan(— ), 



arctan (/3tan(^)) 



(31) 



(32) 



D. 



Q,/3,(T 



COS ( arctan 



(/3tan( — ) 



For a = 1, <; can be obtained from the formula 



2cr 

TT 

Bl, 13,(7, 11 



with 



Bl,i3.cr,ti = + -/9crln(cr) 

TT 



(34) 



(35) 



In the above equations V and W are independent random 
variables, such that V is uniformly distributed in the interval 
(— ^ , ^) an d is exponentially distributed with a unit mean 

oEllll. 

The problem described by Eq. Q and corresponding 
Eq. ( I28t for the stability index a = 2 is a well known 
case of resonant activation resolved in the series of papers 
mmili]. Accordingly, for a = 2 integration in Eq. (|28} 
is performed with respect to a Gaussian measure and standard 
numerical methods can be applied |45| to solve the Langevin 
equation under study. Otherwise, the ordinary forward (or 
backward) in time Fokker Planck Equation |46, 47, 48] for the 
PDF can be derived in a closed form, from which the MFPT 
can be easily calculated. Also, for the Gaussian fluctuations 



theoretical background is provided by the standard thermody- 
namics 1 38] where parameter cr can be related to the intensity 
of the thermal fluctuations imposed on a physical mode x. 

Within these studies, besides the driving Gaussian fluctua- 
tions, two other kinds of a-stable noises have been considered, 
namely the Cauchy noise (a ~ 1, /3 — 0) and the Smirnoff 
noise (a = 0.5, (3 ~ 1) with intensities a = 0.5, a ~ l/\/2 
and (7 = 1. The value of the location parameter /i has been 
arbitrary set to 0. 

The choice of the cr, which scales the distribution width and 
varies the intensity of fluctuations, corresponds, although not 
in a self-transparent way, to the change of system tempera- 
ture. However, the problem of definition of the system tem- 
perature^ is more subtle here than in a presence of a standard 
bath ] lo', 41.'4^.l 50ll5lll because the system is not in the state 
of equilibrium. 

Levy type variables have been generated by use of recipes 
given by Eq. ( I30t and i34l . The trajectory x{t) has been gen- 
erated according to (I28> for various noise realizations. Its 
evaluation proceeded till the time t' when x{t') > 1 for the 
first time. The resulting distribution of first passage times 
(FPT) has been further used to evaluate MFPT. 

The linear potential V± = H±x has been dichotomously 
alternating between different values of H± . Our simulations 
have been obtained for H± = ±8 (cf. Fig. |2] (lower panel)), 
= 0, i7+ = 8 (cf. Fig. 12 (upper panel)) and = 
4, i7+ = 8 (cf. Fig. El(middle panel)). 

Asymptotic lines plotted in Fig.|2]have been calculated nu- 
merically by use of the Monte Carlo method ]52]. As ex- 
pected, typical behavior [29, ,3^ ^3} has been recovered for 
7-> 

(33) MFPT(7 0) = i (MFPT(iJ_) + MFPT(i?+)) , (36) 



and for 7 ^ oo 

MFPT(7-^oo) =MFPTQ(i/_+i/+)^ . (37) 

Fig.|2]presents results of simulations averaged over 10'^ real- 
izations with a time-step At = 10~^. For calculation of the 
asymptotic lines (Eqs. ( I36> and (|^} ) At = 10^^ and averag- 
ing over 10^ realizations has been performed. 

For the barrier switching between H± — ±8, the 
phenomenon of resonant activation is clearly visible cf. 
Fig.|2l(lower panel). By comparison to the results obtained for 
a motion of a par ticle subject to the additive white Gaussian 
noise izSESQlllllIll, flie value of the MFPT is smaller 
and resonant activation is observed at lower frequencies 7. 
Moreover, for H± — ±8, asymptotic values of MFPT esti- 
mated for 7 — * cxD and 7-^0 are higher (lower) than in the 
corresponding white Gaussian noise case I29ll30ll33ll . 

For other cases under consideration, i.e. for H- = 
0, i?+ = 8 and iJ_ = 4, iJ+ = 8, the resonant acti- 
vation has not been observed (cf. Fig. |2] (upper and middle 
panel)). However, values of reported MFPTs for the system 
driven by non-Gaussian Levy noises are always significantly 
smaller than in the case of Gaussian-noise driving. It is caused 
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FIG. 2: The MFPT(7) as a function of barrier fluctuation rate 7 for linear potentials with H+ = 8, H- = (upper panel), — 8, — 4 
(middle panel) H+ — 8, = —8 (lower panel) are plotted along with asymptotic lines. The additive noise is of the Levy-Smimoff type 
(left panel) (a = 0.5, P — 1, fi — 0) and of the Cauchy type (right panel) (q = 1, l3 — 0, fj. — 0), respectively. Different symbols 
correspond to varying intensity of fluctuations: a = .5 (+), a = I/V2 (x) and a = 1 (*). Numerical results have been obtained by use of 
Eq. j28> with the time step dt = 10~^. Noises were generated according to Eqs. <30> and <34> . In order to gain sufficient statistics averaging 
of Eq. l28j over 10'^ realizations has been performed. Asymptotic lines have been evaluated for reference static potential profiles with the time 
step dt — 10^"' and averaged over 10* realizations. Error bars have been estimated by use of standard techniques for the MC data analysis 
1531 . They represent standard deviation from the mean and remain within the symbol size. Lines have been drawn to guide the eye. 



by the fact that the Cauchy and Levy-Smirnoff noises allow 
for higher values of driving fluctuations in the Langevin equa- 
tion Q. 

As it can be inferred from Fig.|2] values of the MFPT for the 
Cauchy noise are higher than the corresponding times for the 
Levy-Smirnoff noise, compare the left and the right panel of 



the Fig. 12] This observation can be explained by the difference 
between both statistics: for the given sets of parameters, the 
Levy-Smirnoff distribution is more "heavily tailed" than the 
Cauchy distribution. 

The procedure applied for numerical integration of ( I28> is 
valid for every allowed value of a. In particular, the case with 
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a = 2 has been investigated along these lines in order to test 
the numerical approach adopted in simulation of Eq. {I). The 
test concluded a perfect agreement of numerically simulated 

results with formerly solved Gaussian cases of RA Eim 

VII. SUMMARY 

We have considered a thermally activated process that oc- 
curs in a system coupled to a non-Gaussian noise source in- 
troduced by a not fully equilibrated thermal bath. Another 
external stochastic process is assumed to be responsible for 
dichotomous fluctuations of the potential barrier which has 
been modeled by the linear function with a varying slope. 

In comparison to the Gaussian case 1 34i|, when the RA phe- 
nomenon has been observed for all barrier setups (H± = 
±8, H+ = 8, = 0, H+ = 8, = 4) analyzed in 
this study, non-Gaussian additive stable noises produce reso- 
nant activation observable only for the H± = ±8 case. Re- 
sulting MFPTs are significantly smaller than those obtained 
for the Gaussian source of fluctuations. The effect is due to 
higher probabilities of the extreme events (large fluctuations) 
allowed by the Levy statistics. Similarly to the Gaussian-bath 
case, a typical asymptotic behavior of the MFPT has been re- 
covered for large and small frequencies of the dichotomous 
noise: for small 7 MFPT tends to the average MFPTs for the 
both barrier configurations, while for large 7 it becomes eq ual 
to the MFPT over the average potential barrier [29, 31, 32]. 

These observations have several implications in relation to 
chemical kinetics in conformationally varying media SEE 
[57'] where flipping barriers separating reactants' and products' 
basins may be due to a dynamic isomerization of the activated 
complex. The effects of fluctuations in force or potential on 
the diffusive process have been also extensively studied in the 
context of motor proteins [58]. Under typical conditions, i.e. 
with a Gaussian additive thermal noise and with a flipping 
barrier height, the RA phenomenon is registered in all those 
realms as a maximal flux of particles or a maximum reaction 
rate. In contrary to this finding, the analogous systems influ- 
enced by the presence of non-Gaussian additive noises exhibit 
RA only for potentials switching between the barrier and the 
well (Fig. |2] lower panel). Otherwise, the kinetics is fairly 
insensitive to the flipping rate of an erecting barrier (Fig. |2j 
upper panel) or becomes slowed down (Fig.|2l middle panel) 
when the frequency of switching between a high and a low 



barrier becomes large. 

In addition, physical interpretation of the Levy noisy term 
in equation Q was discussed in some detail. In general, 
descriptions provided by the FFPE approach and the Tsallis 
statistics are not in a full agreement [5, 36, 37]. In particu- 
lar, the FFPE description of Levy flights in a harmonic poten- 
tial i37l l53il produces different stationary PDF than the one 
predicted by the entropy Eq. ( I19l l. Nevertheless, the Tsallis 
formalism of nonextensive statistical mechanics (or the "su- 
perstatistics" concept, in general [50, 51[) offers an intriguing 
interpretation of non-Gibssian ensembles that can model not 
fully equilibrated baths. Especially, for nonequilbrium sys- 
tems composed of regions that exhibit spatiotemporal fluctu- 
ations of an intensive quantity (like pressure, chemical poten- 
tial, inverse temperature or the energy dissipation rate [50]) 
generalized statistics may emerge in consequence of a statis- 
tical subordination of the intensive variable to those fluctua- 
tions. 

The situation in which resulting probability distributions 
are not Gaussian are frequent in various physical situations. 
The phenomena of special interest are random walks that may 
lead to Levy distributions [5] and hence provide a possible 
realization of Levy noise sources. 

There are several candidates for such a noise origin to be 
taken into account. Obviously, all noise sources need to be 
infinite to allow, with a nonzero probability, infinite fluctu- 
ations. Among possible models of chaotic nonequilibrium 
baths perhaps the most natural is a model of the fluid-like 
bath p?] abruptly perturbed with a local heating to a very high 
temperature. Gradual spreading of energy progresses until a 
new equilibrium state is reached and for some finite periods of 
time, the probing system that experiences local change of the 
bath te mper ature will be governed by non-Boltzman statistics 
i4ll I51I ^5?!. This and similar approaches validate then use 
of the Levy-type statistics for description of not-fully equi- 
librated baths. As known from other similar studies 12811 . the 
Levy-type structure of the noise imposed on a double-well po- 
tential can affect profoundly noise induced jumping between 
metastable states and result in stationary PDF deviating from 
the usual Gibbs distribution. 
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